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ABSTRACT 

We fit the upper main sequence of the Praesepe and Hyades open clusters using stellar models with 
and without rotation. When neglecting rotation, we find that no single isochrone can fit the entire 
upper main sequence at the clusters’ spectroscopic metallicity: more massive stars appear, at high 
significance, to be younger than less massive stars. This discrepancy is consistent with earlier studies, 
but vanishes when including stellar rotation. The entire upper main sequence of both clusters is very 
well-fit by a distribution of 800 Myr-old stars with the spectroscopically measured [Fe/H] = 0.12. 

The increase over the consensus age of —600-650 Myr is due both to the revised Solar metallicity 
(from Zq « 0.02 to Zq « 0.014) and to the lengthening of main sequence lifetimes and increase 
in luminosities with rapid rotation. Our results show that rotation can remove the need for large 
age spreads in intermediate age clusters, and that these clusters may be significantly older than 
is commonly accepted. A Hyades/Praesepe age of —800 Myr would also require a recalibration of 
rotation/activity age indicators. 

Subject headings: 


1. INTRODUCTION 

Stellar models have long been fit to star clusters to 
determine their ages and metallicities; fine grids of evo- 


seem to sugg est an age range of several hundred Myr 
(|Eggen||1998 1, with the main-sequence turnoff giving an 
age of —600-650 Myr for the most massive members 


lYi et al. 2001 

Girardi et al. 

2002 iPietrinferni et al. 


Dotter et al.|2UU8|). l:^'or young clusters, techniques 

indue 

ing fitting models of pre-main sequence contrac- 

tion ( 

Siess et al. 

2000 

iPecaut et al. 2012|) and lithium 

deplei 

ion (|Jettries et a 

.2013 

Hinks &: Jettries|2014 

1 can 


cable to intermediate age clusters, from several hundred 
Myr to several Gyr. Models of stellar evolution are of- 
ten the best dati ng technique available in these cases 
(Soderblom 2010). The ages from stellar modeling are 


(Perryman et al. 1998 Fossati et al. 2008). Other dat- 
ing methods have beeri ap plied to the Hyades, incl uding 
white dwarf cooling tracks ( De Gennaro et al.|2009 ), giv¬ 
ing ages consistent with those from the upper main se¬ 
quence. Ultimately, however, these alternatives rely on 
the same stellar models as ordinary isochrone fitting. 

On a different ground, various authors noticed that the 
presence of extended main sequence turnoffs in 1-2 Gyr- 
old clusters, indic ating an intracluster age spread of sev¬ 
eral hundred Myr (Mackey &: Broby Nielsen 2007 Milone 
et al.|[2009 Li et a!l.||2014 |. 'The ettect of stellar rotation 


t hen used to cali brate secondary age indicators, includ- on the color-magnitilH^iagrams has beeri investigated 


ing those based on stellar rotation and activity (p^ as a possible soMtion to resolve this problem. | Bastiaii 


et al. 1119841 Barnes||2007 IMamajek fc Hillenbrand||2008| 


& de Mink (2009) pointed out that the effects of rapid 


Epstein & Pinsonneault 2014^. 

T he Hyades and Praesepe (the Beehive cluster) are two 
of the best-studied nearby open clusters. As such, they 
currently serv e as benchmarks for stellar ages and stel¬ 
lar modeling (Soderblom 2010). The clusters are both 
part of the Hyades supercluster, and seem to s hare an 
age, as determined both from isochrone fitting (Perry- 
man et al. 1998 Salaris et al. 2004 Foss ati et al. 20081) 
and from stellar rotation (gyrochroiiology, Douglas et al. 
2014 Kovacs et al.|2014 ). The Hyades and Praespe also 
have very similar (possibly identical) chemical co mposi- 


rotation on stellar evolutionary tracks could modify the 
isochrone, leading to a redder and cooler main sequence 
turnoff. However, [Bastian &: de Mink| neglected the ex- 
tension of the main sequence lifetime due to rotation. 


Girardi et al. (2011) pointed out that this change in the 
stellar lifetime i s important, and used the Geneva stellar 


tions (Boesgaard et al.|2013 Taylor fc Joner|2005), hint¬ 


ing at formation in a single molecular cloud or cloud 
complex. 

The age, and even the coevality, of the Hyades and 
Praesepe are not completely settled. Stellar isochrones 
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evolution code (Eggenberger et al. 2008) to get hotter 
and bluer 1.58 Gyr isochrones with a uniform rotation 
rate. These isochrones, like those with core convective 
overshoot but without rotation, did not show a signifi - 
cant extended main sequence turnoff. Yang et al. (2013), 
however, computed their own evolutionary tracks includ- 
ing rotation, and did find extended main sequence turn 
offs in younger clusters. 

While we do not intend to resolve the controversy over 
the extended main sequence problem in this article, those 
efforts hint that rotatio n might be a k ey to resolving the 
age spread observed by Eggen (1998) (Figure 1) in the 
Hyades. Rotation has long been known to have a po¬ 
tentially la rge effect on the evolution ary tracks of stars 
>1.5 Mq (Meynet &: Maeder] 2000). This is distinct 
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from the B5hm- Vetense gap (Bohm-Vi tense 1970) seen 
in the Hyades by de Bruijne et aT 1 2000 ), which occurs at 
the onset of surface convection (arid of efficient magnetic 
braking). 


Recently, Brandt & Huang (2015 ), henceforth BH15 
applied new rotating stellar models (Ekstrom et al.|20T2 


Georgy et al. 2013) to several clusters, including the 
Hyades. I’hey found a best-fit Hyades age of ~800 
Myr from the upper main-sequence turnoff, significantly 
older than the current consensus. We apply the same 
methodology as BH15 to Praesepe as a consistency check. 
We then further investigate whether rotation can resolve 
both the apparent spread in Hyades ages and the incon¬ 
sistency of our older Hyades age with the consensus age 
in the literature. 

We organize the paper as follows. In Section we re¬ 
view our methodology from BH15, referring to that paper 
for a more thorough description. In Sectionj^ we present 
a statistical derivation of the criterion we use to assess 
the consistency of the isochrone-based ages between in¬ 
dividual cluster members. We discuss the selection of 
our cluster samples in Section]^ Sectioncontains our 
results; we conclude with Section]^ 

2. METHODOLOGY 

We apply the Bayesian method used in BH15 to de¬ 
rive the ages of the Hyades and Praesepe open clusters. 
We summarize the method here, and refer the reader to 
that paper for a more thorough description. O ur method 
takes as its input Tycho-2 BtVt magnitudes (Hpg et al. 
2000), parallax, rotational risinj, and a prior probability 
distribution of metallicity, and returns posterior proba¬ 
bility distributions of mass, metallicity, age, and inclina¬ 
tion. In the rest of this paper, we will concern ourselves 
only with the posterior probability distributions of age 
and metallicity. We marginalize over the mass and incli¬ 
nation distributions star-by-star. 


We use t he new rotating isoch rones of Ekstrom et al. 
(|2012|) and |Georgy et al. (|20I3|) as our stellar models, 
i'hese are available at a grid or rotation rates ranging 
from 0 to 95% of breakup, stellar masses from 1.7 to 
15 Mq, and metallicities from Z — 0.002 to Z = 0.014 
([Fe/H] = —0.85 to [Fe/H] = 0). The mass grid is ex¬ 
tremely coarse, so we interpolate these models to higher 
resolution using nonro tating stellar isochr ones. We use 
the PARSEG models (Girardi et al. 2002) for this pur¬ 
pose. We compute the corrections induced by rotation 
as a function of mass, metallicity, and rotation rate, and 
interpolate (and even extrapolate) these coefficients onto 
the much finer PARSEG grid. The validity of this step 
relies on the fact that the rotational correction term is a 
very weak function of the other stellar parameters. 

The rotating stellar models adopt a core overshoot¬ 
ing parameter, the ratio of the convective overshoot to 
the pressure scale height, of 0.1. This value was chosen 
to reproduce the observe d width of the main s equence 
when including rotation (Ekstrom et al. |2012 ). Gore 
overshooting extends the effective size of the convective 
core, allowing more of the star to be burned. Rotation 
achieves a similar effect, but by mixing hydrogen into 
the core rather th a n by simply making the core larger 
(Talon et al.||1997|. Schaller et al. (1992) found that an 
overshoot parameter of U.2 was needed to reproduce the 
observed main sequence width without rotation. This 


degree of overshooting is in mild tensi on with asteroseis- 
mology of two slowly r otating B stars (Papics et al. 


Moravveji et al. 2015). The recent asteroseismic results 
favor an overshooting parameter between 0.1 and 0.2, 
though they find better agreement with exponentially 
decreasing core overshoot than with a step function as 
commonly used in grids of stellar models (including the 
model s we use here) . 


The Georgy et al. (2013) stellar models give the stellar 
equatorial radius, luminosity, and oblateness. We as- 
su me that the star may describe d by a Roche model and 


Espinosa Lara & Rieutord 


(2011) to compute the 


use 

local effective temperature and gravity everywhere on 
the stellar surface. We then use the full specific intensi¬ 
ties of the Gastelli & Kurucz (2004) model atmo spheres 
and the Tycho-2 bandpasses as determined by |Bessell| 
& Murphy (2012) to compute synthetic BtVt photom¬ 
etry as a function of orientati on. The Solar abundances 
have been updated since the Gastelli fc Kurucz (2004) 
models were published, with|Asplund et ai.|(|2U09| End¬ 
ing a Solar metal abundan ce nearly 0.2 dex lower than 
Anders & Grevesse (1989), a nd a Solar iron abun dance 
0.17 dex lower. Because the |Georgy et al.| (|2013| stel¬ 
lar models use the new abundances, we adopt the AT- 
LAS9 atmospheres with [F e/H] = —0.1 relative to the 
Anders & Grevesse (1989) composition for a star. 
This provides an approximate match between the star’s 
bulk and photospheric compositions. It could, however, 
introduce systematic differences in metallicity at a level 
A[Fe/H] ^ 0.1 relative to earlier results, and should be 
kept in mind when comparing different authors’ stellar 
tracks. At an age of ^700 Myr, a decrease of 0.1 dex 
in [Fe/H] can, in some ways, mimic an increase of up to 
^100 Myr in age. 

Finally, we fit the synthetic photometry and predicted 
vsini to the observed photometry, rotation, and paral¬ 
lax. We adopt the a ppropriate priors: masses drawn 
from a Salpeter (19551 initial mass function, random ori¬ 
ent at ions/lmBa uniform prior in volume. We use a prior 
in stellar rotation that closely matches t he observed dis- 
tribution in young, massive field stars (Zorec fc Royer[ 


2012[) . We add a systematic error of 30 kms ^ in v sm i 


and 0.005 mag in Bt and Vr to account for errors in¬ 
duced by our interpolations and finite grid spacing. 

When computing the age and metallicity posterior 
probability distribution of a cluster, we multiply the pos¬ 
terior probability distributions for the individual stars. 
This implicitly assumes the cluster to have a single age 
and composition, i.e., a common origin. We marginalize 
over the other parameters (M, ?;sini, and inclination) 
separately for each star. In the next section, we derive a 
consistency criterion, equivalent to the usual test, to 
assess whether the posterior probability distributions for 
the individual stars really are consistent with a common 
origin. 

3. TESTING THE GONSISTENGY OF POSTERIOR 
PROBABILITY DISTRIBUTIONS 

When determining the posterior probability distribu¬ 
tion in age and composition for a cluster, we multiply the 
probability distributions for the individual stars. Each of 
these distributions is computed using the Bayesian for¬ 
malism detailed in BH15 and summarized in Section 
and marginalized over mass, rotation, and inclination. 
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By multiplying them together, we implicitly assume that 
the distributions are independent in at least one parame¬ 
ter and that they are consistent with a single cluster age 
and metallicity. Our Bayesian formalism, however, does 
not provide a test of these assumptions. 

The assumptions of independence and consistency are 
the same as those used in an ordinary analysis. This 
fact allows us to use the test as a check on the consis¬ 
tency of the probability distributions and, by extension, 
on the ability of the stellar models to reproduce the en¬ 
tire cluster at a single age and composition. We now 
derive an approximate relationship between the product 
of the peaks of the individual distribution functions and 
the peak of their product under the assumptions of in¬ 
dependence and consistency. Unsurprisingly, this ratio 
turns out to be proportional to the likelihood function. 
The y^ test has no effect on our posterior probability 
distributions, but only provides a consistency check on 
the stellar models themselves. 

We begin with the two-dimensional age-metallicity 
probability distributions for each star, the outputs from 
the Bayesian analysis of Section and BH15. We make 
two approximations, well-satisfied in practice, to reduce 
these distributions to one dimension: 


1. The covariance between age r and metallicity 
[Fe/H] is the same for all stars; and 

2. The posterior probability distribution on [Fe/H] is 
entirely determined by our prior. 


In this limit, the constraints provided by the stellar mod¬ 
els are one-dimensional Gaussians extending parallel to 
a single line in T-[Fe/H] space (the covariance matrix 
has one eigenvalue much larger than the other). We de¬ 
fine the variable x to be the linear combination of r and 
[Fe/H] that we are actually constraining, the variable 
that runs in the direction of minimum covariance. We 
obtain the total posterior probability distribution func¬ 
tion by multiplying these one-dimensional Gaussians, 
and then multiplying by the metallicity prior. 

In our case we have A^stars estimates Xi of the actual x, 
where x is a linear combination of age and metallicity (an 
eigenvector of the covariance matrix). For N Gaussian 
measurements of a value Xg, the likelihood £ of x is 


£ oc exp 


(Xi - Xg)^ 


2cr/ 


= exp 


E {Xi - Xg)2 

2a/ 

I ^ 

( 1 ) 


i.e., the exponential of one-half the y^ distribution. Be¬ 
cause the true value xg is not known, the relevant y^ 
distribution has A/tars — 1 degrees of freedom. We can 
then use the standard y^ test to check the consistency 
of the individual posterior probability distributions with 
one another. 


4. SAMPLE SELECTION 


We select our H yades stars from the analysis of [Perry-] 


man et al. (1998), rejecting the blue straggler HIP 20648 
and known spectroscopic binaries, and keeping stars with 


a measured Bt — Vt < 0.3 in the Tycho-2 catalog (Hpg 


et al. 2000). Perryman et al. (1998) rejected several other 


stars that lie above the main sequence. We reject these 


same stars (HIP 20711, HIP 20901, HIP 21670, and HIP 
22565) to enable a direct comparison with earlier results. 
Including all of these targets (apart from the blue strag¬ 
gler HIP 20648) would have almost no effect on our re¬ 
sults. We o btain rotational velocities from the catalog 
compiled by Glebocki & Gnacinski (2005), and a dopt 
the Hipparcos parallaxes and errors for all targets (van 
Leeuwen 2007). Our final Hyades sample includes 14 


stars. 


We s elect our Praesepe sample from the candidates 
listed in Wang et al. (2014) with Bx — Vt < 0.3 (without 
correcting for extinction) and without indications of a 
close binary companion. A fter correcting for E{B — V) = 


0.027 mag (Taylor 


2006), our Praesepe sample has a 


slightly bluer cutoff than our Hyades sample. We exclude 
two stars from this intial sample: the blue straggler 40 
Cnc, and the otherwise unremarkable S Scuti variable BS 
Cnc, which lies significantly below the cluster main se¬ 
quence. This leaves a Praesepe sample of 24 stars. We 
adopt a distance to the cluster of 179±3 pc (Gaspar et al. 
2009| , which we use for each of the individual stars, 'i'his 
is slightly higher than the reported uncertainty of ±2 pc 
(la) to allow for the small, but finite, extent of the clus¬ 
ter core. Most of these stars have risini measurements 


in the catalog of Glebocki & Gnacinski (2005). 

For this analysis, we adopt a metallicity prior of 
[Fe/H] = 0.12 ± 0.04 for both clusters. This value has 
been measured for a l arge sample of Praesepe dwarfs 
(Boesgaard et al. 2013) and is ful ly consistent with th e 
best measurements of the Hyades ([Taylor & Joner|2005[). 
It is slighly higher than the value of [Fe/HJ = 0.10 ±0.05 
adopted by BH15, but is well within the uncertainties. 


5. RESULTS 
5.1. The Age of Praesepe 

Figure [^ shows the results of our analysis for 24 likely 
Praesepe members in the age-metallicity plane. The dot- 
dashed contours enclose 68% and 95% of the posterior 
probability, respectively. The metallicity posterior prob¬ 
ability distribution is almost entirely determined by our 
prior of [Fe/H] = 0.12±0.04. Our marginalized posterior 
probability distribution on age, the red curve at the top 
of Figure gives an age of 790 ± 60 Myr (2a limits). 
This is consistent with the Hyades age given in BH15. 

As for the Hyades, and as discussed in BH15, this old 
age arises from a combination of the increase in main 
sequence of lifetime with rotation, as rotation al mixing 
supplies th e core with additional unburnt fuel (Ekstrom 
et al.|[2M2 ). The blue curve, and lower-left color patch, 
show the posterior probability distribution with the age 
artificially scaled back by this factor. A small additional 
effect arises from the increase in luminosity in the latter 
part of the main sequence. Increasing a star’s luminos¬ 
ity with rotation means that a rotating model overlaps 
a slightly more massive (and shorter-lived) nonrotating 
model in color-magnitude space. 

The best-fit age when fixing the rotation rate in the 
stellar models to zero is much lower: 570 ± 30 Myr. 
This value masks a larger problem with the nonrotat¬ 
ing models: the fitted ages of the more luminous stars 
are inconsistent with those of the less luminous stars, a 
problem that the next subsection will discuss in detail. 
The only way to achieve some measure of consistency 
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Age (Myr) 


Fig. 1.— The age of Praesepe computed by applying the method 
described in BH15 and summarized in Section!^ to 24 early-type 
likely memb ers. We have adopted a metallicity prior of [Fe/H] = 
0.12 it 0.04 lIBoesgaard et al.||2013ll and a reddening E(B — vj = 
0.027 mag (|'i'aylor|2U0b||. Stellar models including rotation indicate 
an age of oiO ± 70 Myr, much older than the consensus age of 
~600 Myr. Much of the increase is due to the lengthening of the 
main-sequence lifetime with rotation. Nonrotating models cannot 
produce a single consistent age, finding a best fit isochrone at a 
metallicity at the extreme end of our prior. 




Fig. 2.— Praesepe posterior probability distributions for subsets 
of the cluster, separated by luminosity. The absolute magnitudes 
assume E{B — V) = 0.027 mag and d = 179 pc. Without rotation, 
the inferred age of the cluster depends on the part of the turnoff re¬ 
gion is used (top panel). The entire turnoff region is not consistent 
with a single age. This inconsistency vanishes when accounting for 
rotation (bottom panel). The entire cluster is consistent with an 
age of 790 it 60 Myr (2 (t), much older than the consensus age. 


is to use a metallicity 2-3(t away from the mean of our 
prior. According to the test, the nonrotating distribu¬ 
tions for the individual stars are inconsistent with one an¬ 
other at 99.99% confidence, while the rotating distribu¬ 
tions disagree at 74% significance, i.e., slightly more than 
the median disagreement expected purely by chance. In 
other words, a single nonrotating isochrone cannot pro¬ 
vide an adequate fit to Praesepe. A rotating isochrone 
with n/Ocrit approximating the observed distribution, on 
the other hand, provides an excellent fit. 

5.2. Hyades and Praesepe Ages by Stellar Luminosity 

With nonrotating stellar models, our analysis for Prae¬ 
sepe fails the test (Section®: the cluster is not consis¬ 
tent with a single isochrone. \Ve now address this prob¬ 
lem by separating our sample into four bins by stellar lu¬ 
minosity (corrected for Ay = 0.084). The bin are chosen 
to have comparable statistical power and to be consistent 
with the Hyades bins we use below. The most luminous 
bin, Vt < 0.9, has three stars, the next bin has four, the 
next six, and the least luminous bin has eleven stars, for 
a total of 24. 

The top panel of Figure clearly shows that the four 
bins of stars are inconsistent with being coeval. The more 
massive, more luminous members require an age of ~550 
Myr (blue curve, close to the canonical age of the cluster), 
while the less luminous members require ages of 600-900 
Myr. The product of the distribution functions produces 
the thick black curve centered near 570 Myr. Simply 
showing this curve masks the strong inconsistency be¬ 
tween the component distributions, which, according to 


the y^ test, have less than a 0.01% chance of being so 
strongly inconsistent by chance. 

The bottom panel of Figure shows that the strong 
discrepancy between the age of the most and least lu¬ 
minous turnoff stars vanishes when accounting for stel¬ 
lar rotation. All stars are consistent with a single age 
and composition. The composition is perfectly consis- 
tent with the spectrosc opic value of [Fe/H] = 0.12 ±0.04 
( Boesgaard et al.|2013 l, but the age, ~800 Myr, is much 
older than the canonical cluster age. This age could be 
somewhat younger when accounting for systematic differ¬ 
ences in the value of the Solar metallicity; at these ages, 
A [Fe/H] = 0.1 can roughly mimic a ~100 Myr difference 
in age. Eliminating core convective overs hoot (currently 
aov = 0.1) from the Georgy et al. (2013) could also pro¬ 
duce a slightly younger age. Asteroseismology of a slowly 
rotating B star, howev er, requires at least th is modest 
degree of overshooting ( Moravveji et al.|[2015 |. 

Figure [^ tells the sarne story as h'igure [^ but for the 
Hyades rather than Praesepe. Because the Hyades is a 
sparser cluster, we use only three bins in absolute Vt 
magnitude. The most luminous bin, Vt < 0.9, has two 
stars, the middle bin has six, and the least luminous bin 
has seven, for a total of 15 stars. These stars extend 
to slightly redder colors than the Praesepe sample (after 
correcting for Praesepe’s extinction). 

As for Praesepe, the Hyades stars are inconsistent 
with a single nonrotating isochrone, though not quite so 
strongly inconsistent. The y^ test gives a 0.4% prob¬ 
ability of at least the observed discrepancy occurring 
by chance. The stars are only consistent with a single 
isochrone at high metallicity, however, [Fe/H] « 0.25. 
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Fig. 3.— Hyades posterior probability distributions for subsets of 
the cluster, sepa rated by luminosity, and only using those stars that 
were also fit by |Perryman et al.| ( 1998||. The top panel shows the 
results without rotation. As tor Fraesepe, the best-fit age depends 
on the stellar luminosity, and the cluster is inconsistent with a 
single age. Also as for Praesepe, the inconsistency vanishes when 
accounting for rotation. 

This is incompatible with the spectroscopic Praesepe and 
Hyades values and the revised Solar composition. Con- 
straining [Fe/H] = 0.12 exactly, a test indicates that 
the stars are fully consistent with a single age when in¬ 
cluding rotation, and inconsistent at more than 99.99% 
probability when excluding rotation. 

5.3. A Consistent, Coeval Picture with Rotation 


More than 15 years ago, Eggen (1998) pointed out that 
the ages of Hyades and Praesepe stars seemed to vary 
systematically with stellar mass. This m ay also be seen 

the 


in Figure 22 of Perryman et al. 
luminous than My 


(1998D: the stars less 
der isochrone than the 


1.8 ht an o 

more massive tur noff stars. This remained true even after 
[Perryman et ah rejected single stars lying above the 
main sequence. The same disagreement between the age 
inferred from the most luminous turnoff stars and the rest 
of the upper main seq uence may be seen in the low er-left 
panel of Figure 17 of David & Hillenbrand (2015). The 
most common solution is to accept the tumott age and 
neglect the modest discrepancy at lower luminosities. 

Stellar rotation, however, enables a fit of the entire 
upper main sequence and turnoff. The previous section 
showed this with stars divided into luminosity bins. Fig¬ 
ure 1^ shows the full color-magnitude diagram, with two 
nonrotating isochrones (600 and 800 Myr, with [Fe/H] = 
0.12) and one 800 Myr, [Fe/H] = 0.12 rotating isochrone 
(yellow and green density plot). The density plot has 
been convolved with the typical observational uncertain¬ 
ties, which are indicated by a cross in the lower-left of the 
image. The units of the density plot are stars per unit Vr 
per unit Bt — Vr- The selection of candidate members is 
described in Section open circles were excluded from 



Fig. 4.— Color-magnitude diagram of the Hyades (blue) and 
Praesepe (red) stars used in our a nal ysis; cluster members excluded 
from the analysis in Sections |5.1| and |5.2| are show n by open circles. 
We ha ve excluded the same Hyades members as [Perryman et al.| 
(|1998|l to enable a direct comparison. Neither cluster is consistent 
with a single nonrotating isochrone, as indicated by the magenta 
curves representing 600 Myr and 800 Myr stellar populations with 
[Fe/H] = 0.12. A rotating 800 Myr isochrone at [Fe/H] = 0.12 and 
with the observed distribution of rotation rates, convolved with 
the typical uncertainty (lower-left cross) and shaded by its pre¬ 
dicted stellar density in the color-magnitude diagram, provides an 
excellent fit to all of the stars. 


the preceding section’s analysis. While increasing the 
metallicity improves the agreement without rotation, we 
find a best-fit [Fe/H] « 0.12 with the rotating models, an 
essentially perfect match to the spectroscopic metallicity. 

A rotating isochrone at a single rotation rate can pro¬ 
vide a fit nearly as good as that from a distribution of 
rotation rates. The necessary rotation rate is very high, 
however, at H/Hcrit 0.7, which is inc ompatible with 


the o bserved distribution of young stars (Zorec & Royer 
20121. This results from the fact that the slower rotators 


have already evolved onto the giant branch, leaving only 
the tail of the distr ibution on the main sequence turnoff 
(Georgy et al. 20141. Adopting an isochrone with a single 
rotation rate would have almost no effect on our inferred 
age, and would modestly degrade the agreement of the 
age distributions of the individual stars. 

The rotating model provides an excellent fit to the ob¬ 
served stellar density in both clusters, confirming their 
common age and composition. The fit is qualitatively 
better than for either of the nonrotating isochrones, and 
accounts for the width of both the upper end of the main 
sequence and of the main sequence turnoff. While the ac¬ 
tual age inferred depends on the value of the convective 
overshooting parameter and on the exact composition of 
the stellar models, our results clearly favor an older age 
for both the Hyades and Praesepe than is currently ac¬ 
cepted. If we adopt an age of ^800 Myr, the apparent 
intracluster age dispersion problem vanishes. 
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6. CONCLUSIONS 

In this paper we have used a Bayesian color-magnitude 
dating technique including rotating stellar models to the 
Praesepe and Hyades open clusters. We have shown that 
the clusters are strongly inconsistent with a single age 
at their spectroscopic composition if we neglect stellar 
rotation. Including stellar rotation, however, makes the 
discrepancy vanish. The Hyades and Praesepe are fully 
consistent with a single episode of star formation and the 
observed distribution of stellar rotation rates. 

However, the age we derive for the Hyades and Prae¬ 
sepe, ^800 Myr, is older than the consensus age of ^600- 
650 Myr. This arises largely from the increase in main 
sequence lifetime with rotation. To a lesser extent, our 
use of models with the updated, les s metal-rich Solar 
composition of Asplund et al. (2009) also increases the 
age. Elim inatiiig the modestc ore convective overshoot¬ 
ing in the Georgy et al. (2013) models could also bring 
our ages closer to the accepted values. However, this 
would likely degrade the agreement of the models with 
observed color-magnitude diagrams, as th e overshoot pa¬ 
rame ter was tuned to match observations (Ekstrom et al. 
2012). Asteroseismic observations also requ ire some de- 
gree of overshooting ( Moravveji e^al.|[2015 ). 

In spite of these caveats, llgurej^shows that the rotat¬ 


ing model reproduces the entire upper main sequence re¬ 
markably well. We have simply adopted and int e rpolated 


the rotating models of Ekstrom et al. (2012); Georgy 


et al. (2013) and fit for an age. We have adopted the 


spectroscopic metallicity and an empirically-motivated 
rotation distribution, adding no additional free parame¬ 
ters to the model. Including rotation removes the need 
for a large spread in ages and suggests that the Hyades 
and Praesepe may be significantly older than is currently 
thought. If the older age is correct, it requires a recali¬ 
bration of secondary age indicators like activity and ro¬ 
tation. 
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